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Abstract Consider N particles moving independently, each one according to a subcritical 
continuous-time Galton- Watson process unless it hits 0, at which time it jumps instanta- 
neously to the position of one of the other particles chosen uniformly at random. The 
resulting dynamics is called Fleming- Viot process. We show that for each N there exists 
a unique invariant measure for the Fleming- Viot process, and that its stationary empirical 
distribution converges, as N goes to infinity, to the minimal quasi-stationary distribution of 
the Galton- Watson process conditioned on non-extinction. 
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1 Introduction 

The concept of quasi- stationarity arises in stochastic modeling of population dynamics. In 
1947, Yaglom [27] considers subcritical Galton- Watson processes conditioned to survive long 
times. He shows that as time is sent to infinity, the conditioned process, started with one 
individual, converges to a law, now called a quasi- stationary distribution. For any Markov 
process, and a subset A of the state space, we denote by \xT t the law of the process at time 
t conditioned on not having hit A up to time t, with initial distribution \i. A probability 
measure on A c is called quasi- stationary distribution if it is a fixed point of T t for any t > 0. 

In 1966, Seneta and Veres- Jones [25] realize that for subcritical Galton- Watson processes, 
there is a one-parameter family of quasi-stationary distributions and show that the Yaglom 
limit distribution has the minimal expected time of extinction among all quasi- stationary 
distributions. This unique minimal quasi-stationary distribution is denoted here u* They 
also show that with an initial distribution fi with finite first moment, fiT t converges to u* s 
as t goes to infinity. 

In 1978, Cavender [13] shows that for Birth and Death chains on the non negative integers 
absorbed at 0, the set of quasi-stationary measures is either empty or is a one parameter 
family. In the latter case, Cavender extends the selection principle of Seneta and Veres- Jones. 
He also shows that the limit of the sequence of quasi-stationary distributions for truncated 
processes on {1, . . . , L} converges to u* as L is sent to infinity. This picture holds for a class 
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of irreducible Markov processes on the non-negative integers with as absorbing state, as 
shown in 1996 by Ferrari, Kesten, Martinez and Picco [TB]. The main idea in [TB] is to think 
of the conditioned process \iT t as a mass transport with refeeding from the absorbing state 
to each of the transient states with a rate proportional to the transient state mass. More 
precisely, denoting N the set of positive integers, the Kolmogorov forward equation satisfied 
by fj,T t (x), for each iGN, reads 

^T t (x)= {q(x,y) + q(x,0) f iT t (y))[fiT t (y)-fiT t (x)}, (1.1) 

where q(x, y) is the jump rate from x to y. The first term in the right hand side represents 
the displacement of mass due to the jumps of the process and the second term represents 
the mass going from each x to and then coming instantaneously to y. 

In 1996, Burdzy, Holyst, Ingerman and March [11] introduced a genetic particle system 
called Fleming- Viot named after models proposed in [18], which can be seen as a particle 
system mimicking the evolution ( II. ip . The particle system can be built from a process with 
absorption Z t called driving process; the position Z t is interpreted as a genetic trait, or 
fitness, of an individual at time t. In the iV-particle Fleming- Viot system, each trait follows 
independent dynamics with the same law as Z t except when one of them hits state 0, a lethal 
trait: at this moment the individual adopts the trait of one of the other individuals chosen 
uniformly at random. Leaving aside the genetic interpretation, the empirical distribution of 
the N particles at positions £ G is defined as a function m(-,£) : N — > [0, 1] by 

1 N 

Wx G N, m(x, - ^ E Mm^}- (1-2) 

i=i 

The generator of the Fleming- Viot process with N particles applied to bounded functions 
f :N N reads 

N oo 

£7(0 = E E id + °) £i m (y> o] [f(C' y ) - /(e)], (1.3) 

%=\ y=l 

where = y, and for j ^ i, C v {j) = an d q{x,y) are the jump rates of the driving 

process. Assume that the driving process has a unique quasi- stationary distribution, called 
z/q S and that the associated iV-particle Fleming- Viot system has an invariant measure X N . 
The main conjecture in [TTIIT2] is that assuming £ has distribution X 1 , the law of the random 
measure m(.,£) converges to the law concentrated on the constant u qs . This was proven for 
diffusion processes on a bounded domain of R d , killed at the boundary [5J Q32 EDI I2S], for 
jump processes under a Doeblin condition [17] and for finite state jump processes pQ. 

The subcritical Galton- Watson process has infinitely many quasi-stationary distributions. 
Our theorem proves that the stationary empirical distribution m(-,£) converges to z/* s , the 
minimal quasi- stationary distribution. This phenomenon is a selection principle. 

Theorem 1.1. Consider a subcritical Galton-Watson process whose offspring law has some 
finite positive exponential moment. Let u* s be the minimal quasi- stationary distribution for 
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the process conditioned on non- extinction. Then, for each N > 1, the associated N '-particle 
Fleming- Viot system is ergodic. Furthermore, if we call its invariant measure X N , then 

VxGN, lim [ \m(x,£) - v*(x)\dX N (£) = 0. (1.4) 

N— >oo J 4 

A simple consequence is propagation of chaos. For any finite set S C N, 

Jim | J] d\ N (0 = Jl^(x). (1.5) 



The strategy for proving Theorem 11.11 is explained in the next section, but there are two 
key steps in the proof. First, we control the position of the rightmost particle. Let 

R(0 ■= . max £(i), 
%e{i,...,N} 

be the position of the rightmost particle of £. Let Q the positions at time t of the N 
Fleming- Viot particles, initially on £. 

Proposition 1.2. There is a time T and positive constants A,cx,C2,C and p, independent 
of N, such that for any £ G 

E(ex V (pR(£)))-ex V (pR(0) < - Cl e" R ^l m>A + Nc 2 e~ CR ^ . (1.6) 

As a consequence, for each N there is a unique invariant measure X N for the N -particle 
Fleming- Viot system. Furthermore, there is a constant k > such that for any N, 
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ex V (pR(0)d\ N (0 < kN. (1.7) 



The second result is that the ratio between the second and the first moment of the 
empirical distribution plays the role of a Lyapunov functional, given that the position of the 
rightmost particle is not too large. For a particle configuration £ define 

: = " " 77^ ■ ( L8 ) 
Recall C N is the Fleming- Viot generator given by (II. 3() . 

Proposition 1.3. There are positive constants v,Ci and Ci independent of N such that 

£>(£) < -^(0 + Ci^ + C 2 . (1.9) 

Propositions 11.21 and 11.31 imply that the expectation of ip under the invariant measure 
is uniformly bounded in N. 

Corollary 1.4. There is a positive constant C such that for all N , 

md\ N (0 < C. (1.10) 
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There are several related works motivated by genetics. Brunet, Derrida, Mueller and 
Munier [HI [TU] introduce a model of evolution of a population with selection. They study the 
genealogy of genetic traits, the empirical measure, and link the evolution of the barycenter 
with F-KPP equation dfU = d^.u — u(l —u) introduced in 1937 by R.A. Fisher to describe the 
evolution of an advantageous gene in a population. These authors also discover an exactly 
soluble model whose genealogy is identical to those predicted by Parisi's theory of mean- 
field spin glasses. Durrett and Remenik [13] establish propagation of chaos for a related 
continuous-space and time model, and then show that the limit of the empirical measure 
is characterized as the solution of a free-boundary integro-differential equation. Berard and 
Gouere [3] establish a conjecture of Brunet and Derrida for the speed of the rightmost 
particle for still a third microscopic model of F-KPP equation introduced in [7J [8] . Maillard 
[2Tj obtains the precise behavior of the empirical measure of an approximation of the same 
model, building on the results of Berestycki, Berestycki and Schweinsberg [3] , which establish 
the genealogy picture described in [TJ [8]. 

We now mention two open problems. The first is to solve the analogous to Theorem II .ll for 
a random walk with a constant drift toward the origin. The second is to obtain propagation 
of chaos directly on the stationary empirical measure, with a bound of order 1/N. 

In the next section, we describe our model, sketch the proof of our main result and 
describe the organization of the paper. 

2 Notation and Strategy 

Let a > and p be a probability distribution on N U {0} such that 



Consider a Galton- Watson process Z t G N U {0} with offspring law p. Each individual lives 
an exponential time of parameter 1, and then gives birth to a random number of children 
with law p. We assume that Galton- Watson is subcritical, that is we ask p to satisfy 



In other words, the drift when Z t = x is — vx < 0. For distinct x,i/6NU {0}, the rates of 
jump are given by 



The Galton- Watson process starting at x is denoted For a distribution \i on N, the law 
of the process starting with /i conditioned on non-absorption until time t is given by 




(2.1) 




(2.2) 



£>-l 




(2.3) 



My) : 



(2.4) 
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where p t (x,y) = P{Zf = y). 

Recall that £f denotes the Fleming- Viot system with generator (11. 3p and initial state £; 
£t(i) denotes the position of the i-th particle at time t. For a real a > define K(a) as the 
subset of distributions on N given by 

Observe that fi G K(a) implies J2xfi(x) < a. 



Proof of Theorem The existence of the unique invariant measure A^ for Fleming- Viot is 
given in Proposition 11.21 

To show (jl.4p we use the invariance of A" and perform the following decomposition. 

m(x,0 - v%{x)\ d\ N {i) = [ E\m(z,$) - u* qs (x)\ dX N ({) 



< \ N (?P >a)+ I E\m(x,g) - v* qs (x)\d\ N (t) 

ip<ot 



< X N ^>a) + / E\m(x,e t )-m(;OT t (x)\dX N (0+ / \m(; £) T t (x) - v* qs (x) \ dX^) 

< X N (ip>a)+ sup \m(;£)T t (x)-v*(x)\ + sup E\m(x,$) - m(;£)T t (x)\, 

(2.6) 

where ip is defined in (11.81) . We bound the three terms of the last line of (12. 6p . 

First term. Corollary 11.41 and Markov inequality imply that there is a constant C > such 
that for any a > 

X N (ip >«)<-. (2.7) 

a 

Second term. Note that "0(0 < at if and only if m(-,£) e K(a). The Yaglom limit 
converges to the minimal quasi-stationary distribution z/* s , uniformly in K(a) as we show 
later in Proposition 17.21 

lim sup \/j,T t (x) - v*(x)\ =0. (2.8) 
Third term. We perform the decomposition 

E\m(x,g)-m{;£)T t (x)\ < E\m(x,g) - Em(x,$)\ + \Em{x, - m{; £)T t (x) \ , (2.9) 
and show that there exist positive constants C\ and Ci such that 

Q e C 2 t 

sup E\m(x,g)-Em(x,g)\ < (2.10) 
££N N y IS 

and c t 

sup \Em{x,g)-m(;Z)T t (x)\ < (2.11) 
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for all N, see Proposition 18.11 later. The issue here is a uniform bound for the correlations of 
the empirical distribution of Fleming- Viot at sites x, y G N at fixed time t. This was carried 
out in pQ. 

To show (11.41) . it suffices to bound the three terms in the bottom line of ( 12. 6p . Choose a 
large and use (I2.7P to make the first term small (uniform in N). Use (I2.8P to choose t large 
to make the second term small. For this fixed time, take iV large and use (I2.9p . (I2.10p and 
( 12. lip to make the third term small. I 

The rest of the paper is organized as follows. In Section [21 we perform the graphical 
construction of Fleming- Viot jointly with a Multitype Branching Markov Chain. In Sec- 
tion H] we obtain large deviation estimates for the Galton- Watson process. In Section [5] we 
obtain large deviation estimates for the rightmost particle of the Fleming- Viot system. In 
Section [6] we study the Lyapunov-like functional and prove Proposition [L3] and Corollary II. 41 
Convergence of the conditional evolution uniformly on K(a) is proved in Section [71 Finally, 
( EDD - dCTD are handled in Proposition Oof Sectional 



3 Embedding Fleming- Viot on a multitype branching 
Markov process. 

In this section we construct a coupling between the Fleming- Viot system and an auxiliary 
multitype branching Markov process (hereafter, the branching process). We call particles 
the Fleming- Viot positions and individuals the branching positions. Each individual has a 
a type in {1, . . . , N} and a position in N. 

When a particle chooses the position of another particle and jumps to it, the process 
builds correlations making difficult to control the position of the rightmost particle. In our 
coupling when a particle jumps, either an individual jumps at the same time or a branching 
occurs at the site where the particle arrives. In this way the particles always stay at sites 
occupied by individuals and the maximum particle position is dominated by the position of 
the rightmost individual (if this is so at time zero). This, in turn, is dominated by the sum 
of the individual positions which we control. 

The coupling relies on the Harris construction of Markov processes: the state of the pro- 
cess at time t is defined as a function of the initial configuration and a family of independent 
Poisson processes in the time interval [0,t]. The coupling holds when the driving process is 
a Markov process with rates {q(x, y), x, y G N U {0}} with being the absorbing state and 
q := sup x q(x, 0) < oo. 

There are two types of jumps of the Fleming- Viot particle i. Those due to the spatial 
evolution at rate q and those due to "jumps to zero and then to the position of particle j 
chosen uniformly at random" at rate q(x,0)/(N — 1). 

Spatial evolution. Each individual has a position in N which evolves independently 
with transition rates (q(x,y),x,y G N) defined by q(x,y) := q(x, y)l{ y ^ so that there are 
no jumps to zero. The spatial evolution of new individuals born at branching times are 
independent and with the same rates q. Under our coupling, each spatial jump performed 
by the 2-particle is also performed by some i-individual. 
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The refeeding and branching. At rate q/ (N — l), each j-individual branches into two new 
individuals, one of type j and one of type i; each new born %- individual takes the position of 
the corresponding j-individual and then evolves independently with rates q. If the 2-particle 
is at x, at rate q(x, 0)/(N — 1) it jumps to the position of the j-particle. Under our coupling, 
each time particle i chooses particle j, each j-individual branches into an i and a j-individual. 
In this way, the i-particle occupies always the site of some ^-individual. 

The branching process has state space 

B:= {CGN^»< N :f;^C(^)<oo} 

i=l x€N 

For i G {1, . . . , N}, x G N, (t(i,x) indicates the number of individuals of type i at site 
x at time t. Let 5fa x ) G 8 be the delta function on (i,x) defined by 8^ x ^(i,x) = 1 and 
8(i,x)(j,y) — for (j,y) 7^ (i,x). The rates corresponding to the (independent) spatial 
evolution of the individuals at x are 

HC, C + % y ) - %*)) = C(h x)q(x, y), i G {1, . . . , N}, x, y G N, 

and those corresponding to the branching of all j-individuals into an individual of type j 
and an individual of type i are 

b {CX + E C(7, x)5 M ) = j^j, /•'./•: {1 V} 

Note that the new born ^-individuals get the spatial position of the corresponding j-individual. 

Harris construction of the branching process Let (JV(i, x,y,k), i G {1, . . . , iV}, x, y G 
N, k G N) be a family of Poisson processes with rates kq(x, y) such that Af(i, x, y, k) C 
J\f(i, x,y,k + 1) for all k; we think a Poisson process as a random subset of R. The process 
x, y, fc) is used to produce a jump of an i-individual from x to y when there are k i- 
individuals at site x. The families (J\f(i,x,y, k), k > 1) are taken independent. Let (J\f(i,j), 
i 7^ j), be a family of independent Poisson processes of rate q/(N — 1), these processes are 
used to branch all j-individuals into an z-individual and a j-individual. The two families are 
taken independent. 

Fix (q = C G B, assume the process is defined until time s > and proceed by recurrence. 

1. Define t(( s ,s) := mi{t > s : t G \Ji !Xi yJ\f(i,x,y,(s(i,x)) U UijjV(i, j)}. 

2. For t G [s, r) define £t = Cs- 

3. If r G A/"(z, x, j/, C s («, x)) then set ( T = ( s + 5^ y) - 5( i>x y 

4. If r G Af(iJ) then set = + ExeN 

The process is then defined until time r. Put s = r and iterate to define for all t > 0. 
Denote the process with initial state £. We leave the reader to prove that Q so defined is 
the branching process, that is, a Markov process with rates b and initial state (. 
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Let |C| := Yli x C(h x ) be the total number of individuals in £. Let 

R(Q : = maxjx : ^ ((x, i) > oj. 

i 

Let Z* be the process on N with rates q and initial position z G N. 
Lemma 3.1. E\C^\ = \(\e^. 
Proof. E\( t \ satisfies the equation 

j t m = ^(EEE^») = ^^(^-1)^101 = ma o.i) 

i x 

with initial condition E\£o\ — I 

Lemma 3.2. Let g : N — > M + be non decreasing. Then 

Eg{R{(<)) < E\$\Eg{Z*®). (3.2) 
Proof. Consider the following partial order on B: 

C^C if and only if ^ ((i, y) < C(h y), forallz,a;. (3.3) 

y>x 

The branching process is attractive: the Harris construction with initial configurations C ^ C 
gives Q -< Q almost surely; we leave the proof to the reader. Let (' := £\ x C(*> x )3(i,R(0) 
be the configuration having the same number of individuals of type i as ( for all i, but all 
are located at r := R((). Hence £ -< (' and 

Eg(R(&) < J29(x)E(Rt,x) < ^gix^Etf&x), (3.4) 

i,x x i 

because g is non- decreasing. Fix i and x and define 

bt(r,x) := Y^E&fax), a t := E\tf\, p t (r,x) := P{Z\ = x). 

i 

Since b t (r,x) and a t pt{r,x) satisfy the same Kolmogorov backwards equations and have the 
same initial condition, the right hand side of ( 13. 4p is the same as the right hand side of (13. 2p . 
This can be seen as an application of the one-to-many lemma, see [2]. I 

Harris construction of Fleming- Viot Let J\f(i,j, x) C J\f(i,j) be the Poisson process ob- 
tained by independently including each r G ftf(i,j) into J\f(i,j,x) with probability q(x,0)/q 
(< 1) by definition of q). The processes (J\f(i,j, x),i,j G {1, . . . , N}, x G N) are independent 
Poisson processes of rate q(x, 0)/(N — 1). 

Fix £ = f e N {1, - ,7V} , assume the process is defined until time s > and proceed 
iteratively from s = as follows. 
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1. Define r(£ s ,s) = inf{t > s : t G U^(z, V> X ) u 6(0)} 

2. For t G [s, r) define & = £ s . 

3. If r G J\f(i,£ s (i),y, 1), then set £ T (i) = y and for i' 7^ i set ^ T (i') = 

4. If r G Af(i,j,£ s (i)), then set £ T (i) = £ s (j) and for i' ^ i set £ r (i') = £ s (i'). 

The process is then defined until time r. Put s = r and iterate to define £ f for all t > 0. 
We leave the reader to prove that £| is a Markov process with generator C N and initial 
configuration £ and the following lemma. 

Lemma 3.3. The Fleming-Viot i-particle coincides with the position of a branching i- 
individual at time t if this happens at time zero for all i. More precisely, 

(0(2, £0(2)) — 1 f or fttti implies (t(i,£,t(i)) > 1 for alii, a.s.. (3.5) 

Corollary 3.4. Assume (o(i,£o(i)) > 1 for alii. Then, 

R(St) < R(Ct), a.s. (3.6) 

4 Galton- Watson estimates 

We show now that for p small enough the functions e p ' belong to the domain of the generator 
of Galton- Watson, that is, the Kolmogorov equations hold for these functions. The total 
number of births of the Galton- Watson process Zf is a random variable H x := x + ^2 t>0 (Z x — 
Zf_) + . Theorem 2 in [22] says that fl2.1l) - fl2.2p are equivalent to the existence of a a' > 
such that 

Eiexpia'H 1 )) < 00. (4.1) 

Clearly a' < a. Let 

F := j/ : N U {0} — > R : J^e^/WI < °o for some p < a'}. (4.2) 

Note that if / G F, then there exist p < a' and C > such that \ f{£)\ < Ce pl , i > 0. For 
/ G F define the Galton- Watson semigroup by 

SJ(x) := E(f(Z*)) < 00, (4.3) 

because < H x for all t > 0. The generator Q of Galton- Watson applied on functions / 
is given by 

00 

Qf(x) := + l)(f(x + £)- f(x)) , x > 0, (4.4) 

t=-\ 

if the right hand side is well defined. 



9 



Lemma 4.1. Under the assumption (12.11) . for f G F ; Qf(x) is well defined and the Kol- 
mogorov equations hold: 

^S t f = QSJ = S t Qf. (4.5) 
at 

Proof. Since |/(x)| < C exp(px) for all ieN, 

\Qf(x)\ < Cxe px (j2p( i + 1 ) ep£ + 1 ))- ( 4 - 6 ) 

This shows the first part of the lemma. Consider / G F and define the local martingale (see 
[21 Section IV-20, pp. 30-37] ) 



/(Z?)-/(x)- f QfiZ:)<ls. 
Jo 



Using gSj), for all s < t 

\M]\ < e p + exp(pH 1 ) +tCH 1 exp(pH 1 ) < CexpipH 1 ), 

with p < p < a'. Hence -Esup s6 [ ^ \M]\ < oo and M] is a martingale by dominated 
convergence. Since for p < a', E exp(pH x ) = (E exp(pH 1 )) x , the same reasoning shows that 
Mf is a martingale and 

Ef(Z?) = f(x) + E J Qf(Z*)ds, 
which is equivalent to (14. 5 p for / G F. I 

The generator of the reflected Galton- Watson process Z t reads 

oo 

Qf(x) := X P( £ + i) 1 ^!} U'( x + *)- f( x )) > ^N, (4.7) 
i=-\ 

if the right hand side is well defined. The reflected process can be thought of as an absorbed 
process regenerated at position 1 each time it gets extinct. Since the absorbed process can 
terminate only when it is at state 1 and jumps to at rate p(0), the number of regenerations 
until time t is dominated by a Poisson random variable Mt of mean tp(0) and 

E(ex V {pZ l t )) < Eexp(pJ2Hl), 

71=1 

where H\ are i.i.d random variables with the same distribution as H l and Mt is independent 
of (H^,n > 1). Hence, 

^(expG^ 1 )) < ex P (tp(0)C(p)). 

Let St be the semigroup of the reflected Galton- Watson process. Using the same reasoning 
as before, we obtain 

Corollary 4.2. Any f G F satisfies the Kolmogorov equations for Q: 

j t S t f = QS t f = S t Qf. (4.8) 
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Large deviations We study Z t , the reflected Galton- Watson process with generator Q 
given by ( 14. 7p . Since p satisfies ( 12. ip . for p < a' < a, 

oo 

T(p) := p(Q) + ^p(£+l)£ 2 e pl < oo. (4.9) 
i=i 

Recall that v is defined in (I2.2p and define /3 as 

(3 = sup{p>0: pT{p) <v}, (4.10) 
which is well defined thanks to the exponential moment of p. 
Lemma 4.3. For any p < min{/3, a'}, and i6N, 

Eexp(pZ?) < e-^ l e px + te p . (4.11) 

Proof. Since p < a' < a, the reflected Galton- Watson generator (14.71) applied to e p ' is well 
defined and gives 

oo 

Q{ef){x) = YsM£ + ±)e px {e pe -l)-p(0)l{*=i}{±-e p ) 

oo 

+ {e p£ - 1 - pi) ) + p(0) l {a:= i } {ef - 1) . 



-i 



£=-i 



Using that for a > 0, e a - (1 + a) < ^e a , 



Q{e"){x) < pxe px (-v + ^T(p))+p(0)l {x=1} e p 

' (4.12) 



; - V 4-e px + e p , 



using p < f3 and /3T(f3) < v. Since p < a', Corollary 14.21 and Gronwall's inequality give 

dun). ■ 

We obtain now a Large Deviation estimate. 

Proposition 4.4. Let T > and 5 > max{l, 4Tp(0)}. Then, there is a constant k, 
independent of x, such that 

j-2 

P(sup (Z: - e~ vs x) >S) < exp ( - 1 — -). (4.13) 

V Q ^r / V rmaxlx.ol/ 



Proof. Set zf = e a; and introduce the process 



e? := - x + v / Zfrfs 



o 



(z*-z x t )+v f \z:- z :)ds. 

Jo 



(4.14) 
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To stop Zf when it crosses 2 max{x, 5} define 

r:=inf|t>0: Zf > 2 max{x, 5}| . (4.15) 

Note that if r < oo, then Z*^ — 2* > 2max{x, 5} — x > 5. Thus, 

{Zf - z$ > 5} C {zf Ar - T > 5} . (4.16) 
For functions gi, g<i : R — > R verifying 



= + u / 92{s)ds, v>0, 



o 



it holds 



Hence, 



Note that 



sup |^i (t) | < - sup |^ 2 (*)| < 5- 



sup |Z* AT - z x tAT \ >5 \ C jsup |e* Ar | >-}. 

t<T J I <<T ^ 









1 sup|e^ Ar | 


>-!} 


-{ 







supef AT . > -1 U {inf e? AT . < — - 
t <T 2 J \t<T - 2 

The treatment of the two terms on the right hand side of the previous formula is similar, and 
we only give the simple argument for the first of them. For p < cr', the following functional 
is a local martingale (see [T5J page 66]). 

M t := exp [pZ x t - px - J" {e~ p Q(e p )) • (4.18) 

Using the bounds of Lemma [4. II we obtain that M. t is in fact a martingale. Observe that 

oo 



P 2 

< -pvx + pp(0) + y(xA(p) + p(0)e p ) : 



(4.19) 



with, 



oc 



A (P) : = 4 E ^ + X ) K - 1 - A*) > 0. ( 4 -20) 



We have already seen that A(p) < T(p). Then, we bound the martingale M. t as follows. 

^2 /•< ^2 



M t > exp (p(Zf - z ) - ( - p„ + ^-A(p)) jf Z> - pp(0)t - P -te" 
> exp (pe* - pp(0)t - ^r(p) jf Z> - ^te" 



(4.21) 
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By stopping the process at r, and using that 5 > 1, we obtain for t < T 

exp (pe x tAT ) < Mat exp (pp(O)T + p 2 max{x, <f}7T(p)) . (4.22) 
Using (glSD, (Q7|) and (Q2|l . and the bound p(0)T < 5/4, we obtain for any p > 

P(sup - z*) > 5) < p( supM sAr > exp (£- - p 2 max{x,5}Tr(p)) N ) 

s<T ^ <s<T 4 / 

P 5 (4 ' 23) 

< expf — — + p 2 max{x, 5}Tsup r(p) ) , 

by Doob's martingale inequality and for p < (3. Optimize over < p < (5 (recalling that 
16T/3T(/3) > 1), and choose p* 

< (3. (4.24) 



16T max{x, 6} T(/3) 

The result follows now from ( I4.23P and (I4.24p . I 

5 Bounds for the rightmost Fleming- Viot-particle 

In this section, we bound small exponential moments of the rightmost Fleming- Viot-particle. 
We first define a threshold A, such that with very small probability, the rightmost particle's 
position does not decrease when it is initially larger than A. Define 

^4( 1 - exp (-iii))) €(0 ' 1) ' 

Choose 

min{/3, a', 7Kp(0)} 
Po := ~ A 



where k is the constant given by Proposition 14.41 Define 



a=**<°>>i. (s.i; 

Po 



Define the time and the error 5 entering in the large deviation estimate of Proposition 14.41 
as follows. For an arbitrary initial condition £, 

T:= 4p70)' S: = max {^^J~}^ ( 5 - 2 ) 

recall here that R(£) = m&Xi< N and set Vl(£) = exp(pmin(.R(£), L)) for L > A which 
will be taken to infinity later. We use the notation [-F(£t)]o := F{£.t) — F(£,o)- 
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Proof of Proposition We use the construction in Section [3] to couple the Fleming- Viot 
process £f and the branching process Q with ( = £)t so that ((i,l;(i)) = 1 for all i. 
Then, by (13. 6 \ R(£t) < R(Ct) an d h is sufficient to prove an inequality like (11. 6p for R((t)- 
Notice that for the initial configurations £ and £, -R(£) = -R(C)- We drop the superscripts £ 
and £ in the remainder of this proof. 

Define the event 

G = Q{& T) := {R(Ct) - e~ vT R(0 < 5} , (5.3) 
and for a positive real c, we define the set 

K c :={t:R(£)<c}. 

On K% 5 = R/A < R, and on K% n 0, 

«(Ct) < (1 + ^)^) < (1"7)«(0- (5-4) 

Hence, 

lj^no [Vt(C*)]J < Vt(0 (e-^ - 1) lK>nK L no < ~Vt(0 (l - e"™) l^nK.ng. (5-5) 
Since A > 1, on AT A n £, #(Ct) < ^e"^ + 1 < 2 A so that 

Thus 

[V L (Ct)fo < -(l-e-^)e^)lKcnK,n 6 + e^l^no + K (C) ]o V 

< _(i_ e -7P^ yL(e)1 ^ n ^ + e 2 Ml ^ + 2e ^)l C c, 

where we used that 



(5.6) 



Choose p := min(p , 47^42) an d observe that by Lemma [372]. 

E [e 2 » R ^] < E\(t\ E [exp(2pZ? e) )] < Ne^ T ( e -W JpW) + Te 2p ) . 

by Lemma 13.11 for the bound of the first factor and Lemma 14.31 for the bound of the second 
factor. Also, Lemma [4.41 implies 

P{Q C ) < £|Cr|Pf sup (ZjW - e~ vs R(0) > s) < Ne^ 7 \e~^l m<A + e~^l m>A ). 

v s<T ' 

Taking expectation on (15. 6p we bound the last term as follows. For constants C\,Ci,C\ y and 
C 2 



E[e"^l m ] < (P(g c )E[e 2pR ^]) 1/2 



< C 1 Nl KA + C 2 Nexp(- KR ^' 



1/2 



k a (5.7) 
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Gathering ( 15. 7p and ( 15. 6 p we obtain, for any L > A, 

EV L (&)-V L (0 < -c 1 V L (0l L >R(0>A + C 1 Nl m ^ A + C 2 Ne-"~ C2R ^ 

< -cMOlLyRi^A + CNe-^ 

which completes the first part of the proof, inequality (11.61) . as one takes L to infinity in 

(EHD. 

For the second part, take C > and observe that the set of £ such that the right hand 
side of (15. 8 p is larger than — C is finite. Foster's criteria, [23| Theorems 8.6 and 8.13] implies 
that both the chain (£o ; £t,£2T, • ■ •) and the process ^ are ergodic with the same invariant 
measure that we call A". 

Now, consider again (15. 8p for a fixed L. Note that Vl is bounded, so that by integrating 
(15. 8p with this invariant measure, and then taking L to infinity, we obtain (jl.7p . ■ 



6 The empirical moments of Fleming- Viot 

In this section we prove Corollary 11.41 Introduce the occupation numbers i|:Nx N w — > N 
defined as 

N 

8=1 

for which we often drop the coordinate £. Notice that m(x, £) = rj(x,£)/N. 
For any integer k, define the fc-th moment of the N particles' positions as 

N oo 
i=l x=l 

As there are only N particles, M k is well defined. Instead of working with the barycenter 
Mi/N, we consider ip := M<ijM\. Note the inequalities 

i < ^ < m < mi (6.i) 

The function ip is not compactly supported (nor bounded). Even though C N ip is well defined, 
we need to use later that J C N ipd\ N = 0. We do so by approximating ip by a compactly 
supported function tp L for which we have 

C N ij L d\ N = 0, and lim £> L = £> pointwise. (6.2) 

L— >oo 

We approximate the unbounded test function ip by the following one 
M L (f) N L 

ri0 = M^y with ^(0 = E™^ fc ( < )^ fc ) = E a;fc ^,0+^E^0- (6-3) 

1 i=l x=l x>L 
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As N is fixed, = L k N — J2^. =1 (L k — x k )r](x), and has compact support. It is easy, and 
we omit the proof, to see that there exist a positive constant C such that 

|£> - £> L | < |£>| + |£> L | <Ci><CR, (6.4) 

where we recall that R(£) = maxj£(i). We have established in Proposition 11.21 that R(£) is 
integrable with respect to A , so that (16. 2p implies that 

£ N ^d\ N = 0. (6.5) 

The main result of this section is the following. 

Lemma 6.1. There are positive constants C±, C2 such that for any integer N large enough, 

^2 f d2j\n 



^d\ N < + ^ J R d\ N . (6.6) 

Proof of Lemma 16.11 We decompose the generator (II. 3p into two generators, one govern- 
ing the refeed part and the other the spatial evolution of the particles: C N = C^nit + ^refeed' 
which applied to functions depending on £ only through 77 (•,£), read 

Cfeed=P(OMl)f;^ T K^-l), With A£{ V )(y) = V**, ( g <7) 

~l N - 1 [vw ±1 y = x, 

00 00 00 

Cift = E^M )^^-! - 1) + Y^xnM + !)(^^ " !)• (6-8) 

x=2 x=l i=l 

It is convenient to introduce a boundary term 

B = —r](l)p(0)(Ai A£ — 1) and call £q = £^ rift - B, (6.9) 
which applied on -0 yield 

, s , s / M 2 - Mi \ , 

x=l i=— 1 x 7 

^ , . , ^f2xM 1 -M 2 + iM 1 \ 

= E^)[E^ + 1) ^ Ml (M 1+ J 

, . M 2 - Mi / x x - . Mi \ M 2 



(6.11; 



-x ■ - i = l 



00 



< + p(0) ^ + J2p(i + l)i 2 < -W + Co 
1 %=\ 
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for some positive constant Cq. Finally, for the jump term 

Aefeed^ = p(0)^(l) £ J^J { ^ + x _ , ~ ^ ) 

X Z (6.12) 
- p( ° M1) ^iV-l M^M-l) X l + ^' 

If we set A(x) = 1/(1 + x) - (1 - x), for x G [0, 1], then 

A(x) = , and < A(x) < x 2 . (6.13) 

We apply (I6.13P to expand the last term in (I6.12p . with x/(M\ — 1) < 1 for x < -R(£); and 
obtain 

L N ib- v(0)n(l) V 7/(z) M(x 2 -1)-M 2 (x-1) / x x A 

^efecd^-P(0Ml)^ iv _ 1 M^M.-l) *[} A/ 1 -l +A( Af 1 -lV' 



(6.14) 

Note that 

oo oo 

^2 {Mtx 2 - M 2 x) = 0, and ^ r/(x) (M^ 2 - M 2 x) (-x) = -M 3 M 1 + (M 2 ) 2 . 

x=l x=l 

Also, 



^ iv - 1 v ; V Mi - ly V Mi - 1/ 



:(M 2 -Ml) 



Mi - 1 y N — 1 
1 

(JV-lXMx-l) 

M 2 - Mi 



(M 2 - Mi) 



(iV-l)(Mi-l)' 
Thus 

„ w / x 77(1) M 3 Mi - (M 2 ) 2 . . . . M 2 - Mi 

where 

_ p(0)r / (l)(M 2 - Mi) ^ ^) Mi(x 2 - 1) - M 2 (x - 1) x 

KeSt " (iV-l)(Mi-l) + ™ 1] ^JV-1 Mi(M!-l) XA( M 1 -1 ) - 

(6.15) 

Using that M 2 - M 1 > 0, 
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Thus, we reach that for C independent of N and L, 

C N W < -vip + 24p(0)^ + C . (6.18) 

We now integrate (16.181) with respect to the invariant measure, and use that J C N ipd\ N = 
to obtain for constants Ci, and C<± (independent of N) 



I 



f R 2 d\ N 

^d\ N <C 1 + C 2 J ■ (6.19) 



7 Uniform convergence to the Yaglom limit 

Define the generating function of a distribution /i on N by 

G(wz) := ^2^(x)z x , zeR,\z\<l. (7.1) 



In this section we show a uniform convergence for \x € K(a) of the generating functions of 
fiT t to the generating function of the qsd v. We invoke a key result of Yaglom [27] . The 
continuous time version can be found in Zolotarev [28] . 

Lemma 7.1. [Yaglom 1947, Zolotarev 1957]. There is a probability measure v such that 

lim G{8 1 T t ] z) = G{v] z), (7.2) 

t— >oo 

and the generating function of v is given by 

G („ ; ,) = l-exp(-„£ E — ^-), , eP ,l). (7.3) 

The measure v is in fact u* s , the minimal QSD. We do not use the explicit expression 
( 17. 3 p of the generating function of v\ we only use (17.20 . Recall that [iT t is the law of Z t with 
initial distribution \x conditioned on survival until t and that K(a) is defined in (12.51) . The 
next result says that the Yaglom limit holds uniformly for all initial measures in K(a). 

Proposition 7.2. For any a > 

lim sup \G(fjTt; z) — G(u*; z)\ = 0. (7.4) 

As a consequence, for each iGN, 

lim sup \iiT t (x) — v* (x)\ — 0. (7.5) 



t— >oo 



A' (a) 
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Proof of Proposition 7J}. Recall that St is the semigroup of the Galton- Watson process and 
observe that for any £ G N, G(5iS t ] z) = G l (5iS t ] z). We set, for simplicity, 

g(z) := l-G(S 1 S t ;z) G [0,1], 

for z G [0, 1]. The following inequalities are useful. For z G [0, 1], 

l-£g(z) < (l-g(z)Y < 1 - £g(z) + fg 2 (z). (7.6) 

The generating function of fiT t reads (the sums run on I G N) 

G{\iS t \ z) - G(iiS t ; 0) 



G(nT t ; z) 



1 - G(fiS t ; 0) 
Z e M(G(5 e S t ;z)-G(5 i S t ;0)) 

E,/iW(l-G(^;0)) (7-7) 

E^W((i-y(^- (i-»(o))') 



Also, 



1 - G{5 x T t] z) 



e^w(i-(i-^(om 



l-G(5 1 S t ;0) g(0Y 

We now produce upper and lower bounds for G(fjT t ; z) — G(z/* s ; z). We start with the upper 
bound. Using first (I7.7|) and then 07.6p . 

J2 if i(£) ((1 - g(z)Y - 1 + (1 - (1 - g(0)Y)(l - G« s ; z))) 
GW;z)-G K ; Z ) = E ^)(l-(l-^) 

£< W ( - »W + + 0(O)(1 - G« s ; «))) 



< 



EM£){g(o)-tg 2 (o) / 
EMQ(0- - G K; *)) - fi) + E^VWf^W 



< 



£/4*0(i-4Kt») 



1 _ M2 ^qfQ') 

1 Mi( M )yw 



(7.* 



where M fc (/x) := E^VW, fceN. Thus, 



^ T v r( . , . \G(6 1 T t -z)-G(v* s] z)\ + (l-G(5 1 T t -z))g(z)a 
sup G{jjT t ;z)-G(y ;z) < ^ — . (7.9) 
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Now, for the lower bound, we use similar arguments to reach 

Ee W) ( - g{z) + »(0)(1 - G(v* s ; z)) - V(0)(1 - G« s ; z))) 



G(fxT t ;z)-G(v*-z) > 



> G(6 1 T t - z) - G(»* qa ; z) - v^(0)(l - G« s ; z)). 



MM 

M 1 (/i) i 



(7.10) 



Thus, 



inf G{pT*z)-G{v^z) > -\G(6 1 T t] z) - G(i£; z)\ - a»(0)(l - G(i£;z)). (7.11) 

Since g(z) goes to as the implicit i goes to infinity, both (17. 9p and (17.111) go to 0. This 
proves (I7.4p . The proof of (12. 8p follows from (17.41) and Lemma [7.31 below on convergence of 
probability measures. I 



Lemma 7.3. Let {fi^,n £ N, 7 £ T} be a family of probability measures 
each z £ [0, 1] we have 

lim sup \G(f/l, z) — G(v, z)\ =0. 

Then, for each xeNwe have 

lim sup \[J>n( x ) ~ v i x )\ = 0- 

Proof. Let / = We consider the one-point compactification of N, which we denote 

N = NU {00} and extend /: N — > K by /(oo) = 0. Since / is continuous function on N, the 
Stone- Weierstrass approximation theorem yields a function h, which is a linear combination 
of functions of the form {y 1— y a y , < a < 1} (finite linear combinations of these functions 
form an algebra that separates points and contains the constants), such that for any e > 0, 

sup^N \f(y) - %)l < £ - Then 

sup\/il(x) - u(x)\ = savlulf -vf\ < sup - + sup \^ n h - vh\ + \vh-vf\. 

7 777 

The first and the third term on the r.h.s. are smaller than e while the second one goes to 
zero as n goes to infinity by assumption. ■ 



Assume that for 
(7.12) 



8 Closeness of the two semi-groups 

In this section we show how propagation of chaos implies the closeness of Em(x,^) and 
m(-,£)T t uniformly in £ £ A^. The arguments are similar to those used in [I7J Q]. The key 
is a control of the correlations that we state below. For a signed measure /1 in N we will 
need to work with the £ 2 norm given by ||/x|| 2 = ^2 x& ^(fi(x)) 2 ■ 
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Proposition 8.1. There exist constants c and C such that, 



Ce ct 

sup \\E[m{x,g)]-m(;t)T t \\ < — . (8.1) 

?eN AT 1\ 



As a consequence, 



Ce ct 

sup \E[m{x,$)]-m(;Z)T t {x)\ < — -. xeN. (8.2) 

Furthermore 

2 ' 

sun F!\m,(T..£.} — m,(- P\TA < - 

N 

Proposition 8.2 (Proposition 2 of pQ). For eac/i £ > 0, and any x,y 6 N 

sup |F[m(x,^)m(y,^)]-FH^d)]^H^ ; ^)]| < • (8-4) 



sup E[m(x,$) -m(-,0T t ] 2 < xGN. (8.3) 



The paper [T] proves this proposition for processes with bounded rates, but the extension 
to our case is straightforward. 

Proof of Proposition ^. 1\ Fix £ G and introduce the simplifying notations 

u(t, x) := Em(x, £|) and t>(t, x) := m(-, £)T t (x). (8.5) 
Define x) = u(t, x) — v(t, x). We want to show that for any t > 0, 

liiwir < |ii^)ir+ 4p( y o)f - (8.6) 

Recall the definition (12. 3p of the rates q and the evolution equations satisfied by v(t,x) and 
u(t, x): 



■^v(t,x)= ^ q(z,x)v(t,z) - (^2q(x,z) > jv(t,x) + p(0)v(t, l)v(t,x), (8.7) 

zj^x,z>0 z^=x 



—u(t,x)= Yl q{z,x)u(t,z)-(^q(x,z)y(t,x)+p(0)u(t,l)u(t,x)+W(^,t,x). (8.8) 

zj^x, z>0 z^x 

Here, 



W{i- 1, x) = p(0) (-JL^Elmfa £)m(l, £)] - £[m(l, £[m(x, £)]) . (8.9) 
Proposition 18.21 implies that 

su V \W&t,x)\<^- . (8.10) 

€ 7V 
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Observe two simple facts. First, set D = {(x, z) : x > 1, z > 1, x ^ z}, and for any function 
f :N^R 

{q(x,z)+q(z,x))f 2 (x)-2 £ q(x, z)f(x)f(z) = ]T - f(z)) 2 . 

(x,z)eD (x,z)eD (x,z)eD 

(8.11) 

The second observation is specific to our rates. For x > 

$>(2,x)<^g(x, 2 )+p(0). (8.12) 

ZJ^X Zy^X 

Observation ( 18. lip is obvious and we omit its proof. Observation (I8.12p is done in details. 
^^q(z,x) = zp{x — z + 1) = x p(x — z + 1) + (z — x)p(x — z + 1) 

Z^X Z>0,ZJ^X z>0,zj^x z>0,z^x 

= x(p(0) + p{\) + ■ • • + p(x + 1)) + (p(0) - p(2) xp(x + 1)) 

< xj^p(z') +p(0) = 2) +p(0). 

(8.13) 

Now, we have 
d 

^2S(t,x)—5(t,x) = ^2 (q(z,x)5(t,x)5(t,z) - q(x,z)5 2 (t,x)) 

x>0 (x,g)eD 

+ p(0) ^ (u(t, x)u(t, 1) - x)v(t, l))6(t, x) + J2 x)W{& t, x). 

x>0 x>0 

(8.14) 

Let us deal with each term of the right hand side of (18.141) . For the first term we use (18. lip 
and fl8TT2|) . 



(x,z)eD 

< q(z,x)6(t,x)6(t,z) -^^(^g(x,z) + ^g(2;,:r) -p(0))5 2 (t,x 

(x,z)&D x>0 z^x z^x 

<~ E ^,x)(5(t,x)-5(t,,)) 2 + ^ii<5(t)ii 2 

(a:,z)eD 

/ P(0)nx/.Mia 



2 

(8.15) 



To deal with the second term, first note that 



sup 1 < = iwoii- 



z>0 



x>0 
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Then, 

^2(u(t,x)u(t,l) -v(t,x)v(t,l))5(t,x) < (S(t,x)u(t, 1) + v(t,x)5(t,l))6(t,x) 

*>o x->o (gl6) 

< ^ x) + m, i) i sup \s(t, x) | x ) < 2 \\m f- 

x>0 x>0 x>0 

For the last term, we have 

\J2Ht,x)W(^t,x)\<sup\W(e,t,x)\xJ2\Ht,x)\ <2s* V \W(Z;t,x)\. (8.17) 

_ „ x>0 . _ x>0 

x>0 x>0 

Thus, we obtain (18.61) . Gronwall's inequality allows to conclude. I 
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